Impurity effects in a two— dimensional system with Dirac spectrum 
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It is demonstrated that in a two-band 2D system the resonance state is manifested close to the 
energy of the Dirac point in the electron spectrum for the sufficiently large impurity perturbation. 
With increasing the impurity concentration, the electron spectrum undergoes the rearrangement, 
which is characterized by the opening of the broad quasi-gap in the vicinity of the nodal point. If 
the critical concentration for the spectrum rearrangement is not reached, the domain of localized 
states remains exponentially small compared to the bandwidth. 

PACS numbers: 71.23.An, 71.30.+h 
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The effect of impurities on quasi-particle spectra in 
disordered systems is qualitatively determined by the ra- 
tio of the dimensionality of the disordered system to the 
exponent in the respective dispersion relation. The dis- 
ordered system exhibits low-dimensional behavior when 
this ratio is less than unity. In general, impurity effects 
are more pronounced in low-dimensional systems. Ma- 
terials like graphene are certainly 2D objects^. How- 
ever, electrons in graphene feature the linear dispersion 
close to the Fermi level. A number of experiments ev- 
idently demonstrate that graphene is highly tolerant to 
impurity induced perturbations. This fact can be at- 
tributed to the increased effective dimensionality of the 
electron subsystem in graphene. With respect to the or- 
dinary quadratic dispersion graphene could be regarded 
as a four-dimensional system. This high effective dimen- 
sionality should be beneficial for the reduction of local- 
ization effects that occur due to impurities, which are 
inevitably present (or intentionally introduced) in corre- 
sponding materials. 

The importance of impurity effects for the physics 
of graphene had been frequently emphasized. Notwith- 
standing, the effect of disorder were studied only in both 
weak scattering^ and unitary^ limits, or for a kind of 
interpolation between such extreme cases&. When im- 
purity states of single defects are located in the vicinity 
of the van Hove singularities of the host system, an in- 
crease in the impurity concentration yields a substantial 
spectrum rearrangement (SR), albeit the relative impu- 
rity concentration remains quite low££. This transition 
between two qualitatively different regimes of impurity 
scattering takes place only for a finite magnitude of the 
single-impurity perturbation. The type of the state that 
is produced by the single impurity is usually reflected in 
the passage of the SR. Below we are attempting to ex- 
amine a possibility for impurity states to appear close to 
the Dirac point of the electron spectrum in a 2D system 
with linear dispersion for an arbitrary strength (unitary 
limit including) of the single-impurity perturbation, and 
to outline a scenario of the SR with varying the impu- 
rity concentration. Similar issues have been raised up in 
Ref. but the problem have not been solved correctly. 



In order to model a system with the Dirac spectrum, 
one can choose the host tight-binding Hamiltonian in the 
most basic formifl, 

h = E[/( fe H( fe ) c 2( fc ) + r(fc)4(*0ci(fc)L (i) 

k 

where c^(fe) and c a (k)) are creation and annihilation op- 
erators on two sublattices, and k is a 2D wave vector. 
Since only the close vicinity of the nodal point will be of 
concern, it is sufficient to put f(k) = ta(k x +ik y ), t > 0, 
where t is the hopping parameter, and a is the lattice con- 
stant. Then, the dispersion relation, e(fc) = ±tak, does 
possess a Dirac point at the zero energy, which separates 
two bands that are touching each other. 

We also assume that our system can be reasonably well 
described as a substitutional binary alloy with a diago- 
nal disorder (a so-called Lifshits model). It is supposed 
that impurities are distributed absolutely at random on 
both sublattices, so that on-site potentials can take one 
of two values, say Vl and 0, with probabilities c and 1 — c, 
respectively. The full Hamiltonian of the disordered sys- 
tem is then represented by the sum of the translationally 
invariant host part Q and the perturbation, 
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where < a,p > ranges over those sites on the lattice that 
are occupied by impurities. 

Let only the zeroth site on one of sublattices be occu- 
pied by an impurity. Then, the diagonal element of the 
Green's function (GF) G — (e — H)^ 1 on this site 



G = go/0- - v L9o), 



(3) 



where go is the diagonal element of the GF in the host, 
g = (e — Ho) ■ Site-diagonal elements go are equal on 
both sublattices and can be easily obtained by approxi- 
mating the Brillouin zone with a circle, 
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FIG. 1: LDOS at the impurity site for v = -10, -5, -2.5 is 
shown by solid, dashed and dot-dashed curves, respectively. 
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It is convenient to choose the energy unit in such a way 
that the bandwidth = 1. Thus, for J3J one obtains 



.g = eln(e 2 /(l-e 2 ))^|e| 



(5) 



The local density of states (LDOS) at the impurity site 
(see Fig.P) is given by the imaginary part of the diagonal 
element of the GF 
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(1-veln [e 2 /(l-e 2 )]) 2 + (z;7re) ; 



r, (6) 



where v is the dimensionless single-impurity perturba- 
tion. For the sufficiently large |u|, a prominent peak is 
manifested in the LDOS l@ close to the Dirac point in 
the spectrum, indicating the presence of the resonance 
state. Its energy e r is defined by the Lifshits equation 



1 » 2ve r In |e r 



(7) 



It should be emphasized that for the attractive impurity 
potential v < the energy e r is located above the nodal 
point (e r > 0), and, vice versa, it is located below this 
point (e r < 0) for the repulsive impurity potential v > 0. 
In contract to 3D systems, the resonance state is accom- 
panied by the deep local level outside of both adjacent 
bands. Thus, the total number of states near the nodal 
point is gradually diminishing with increasing \v\. 

When the resonance peak is relatively narrow, the de- 
nominator in Eq. 10 can be expanded about e r , 
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The resonance state is well-defined when the effective 
damping JSJ is much less than its separation from the 
closest van Hove singularity, 



r/|e r | W7r/[2|ln|e r | + l|]<l. 



(9) 



The inequality is satisfied only when the resonance en- 
ergy e r is located fairly close to the Dirac point and is 
strengthening with decreasing |e r |. Thus, the resonance 
presence in the unitary limit^ is justified. It is worth 
mentioning that well-defined resonances can not appear 
in the vicinity of the band edge in a single-band 2D or 
3D system within the Lifshits model J2J- 

It is not difficult to calculate also the change Ap in 
the total DOS in the system that is caused by the single 
impurity center—, 
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It can be verified that bare states are redistributed within 
the bands. For the case v < 0, states are removed from 
the domain of the continuous spectrum in the lower band 
to the split off local level, and in the upper band states 
are pushed towards e r . However, there is a notable neg- 
ative dip in Ap at the nodal point, where the host DOS 
is zero (see Eq. (0). Therefore, the close vicinity of the 
nodal point can not be described properly by the direct 
expansion in the impurity concentration even at the neg- 
ligibly small concentration of impurities. 

Commonly, renormalized methods, such as the cohe- 
rent-potential approximation (CPA), are the most effec- 
tive inside the continuous spectrum. The one-electron 
GF of the disordered system can be expressed by the 
corresponding self-energy S(fe). Since the translational 
invariance is restored by the configurational averaging 
over impurity distributions, 



G(fc)" 1 = g(k)- X -S(fe), 



(11) 



where operators G(k), g(k), and fi(fe) are acting in the 
sublattice space. For the model system under consider- 
ation, the self-energy within the CPA is site-diagonal 
and identical on both sublattices. According to the con- 
ventional procedure, it should be determined in a self- 
consistent manner from the equation, 



cv/(l -(v- cr)g Q (e - a)). 



(12) 



In the effective medium constructed by the CPA, the 
self-energy can be expanded into the series in impurity 
clusters^, 

£^(fc) = 8 al} ci + af{k) + ... 1 (13) 

where ^(fe) represents the contribution from pair dia- 
grams, 

- a0 (k) = 6 af3 
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In Eq. (|14f> I and indices a, f3, and 7 enumerate lattice 
cells and sub-lattices, respectively, 
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the single-site T-matrix is denoted by 

T m = [V m - f)/[l - ( v rn - CT).9o(e - O")], 



(16) 



while indices m and n enumerate atom types (impurity 
or host), so that £ m attains values c or 1 — c depending 
on the value of these indices, and variable v m is v or 0, 
respectively. 

The relative magnitude of contributions from scatter- 
ings on impurity clusters is increasing as approaching any 
van Hove singularity in the spectrum, so that the CPA 
becomes unreliable in their vicinity. The necessity to im- 
plement a relevant applicability criterion for the CPA and 
other approximate methods based on the partial sum- 
mation of the series for the GF have been overlooked in 
some recent articles devoted to the impurity effects in 
graphene 5 . The analysis of the series expansion for S(fc) 
shows that the series does have a small parameter, 
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Cluster diagrams can be omitted on |i?(e)| ^ 1/2. Inside 
the energy domains, where this inequality holds, only the 
first term can be retained in the series and the resulting 
approximate expression for the self-energy does not de- 
pend on fc. If the relative impurity concentration is kept 
low, multiple-occupancy corrections that are included in 
the derivation of the CPA can be neglected too, so it is 
reduced to the so-called method of modified propagator, 



cv/(l ~ vg (e - a)). 



(18) 



Since our interest is restricted to the narrow vicinity of 
the nodal point in the spectrum, it is possible to make 
an obvious approximation for the diagonal element of the 
host GF, 



go sa 2eln|e| -wr|e| , |e| < 1. 



(19) 



By making a substitution e — a = Hexp(itp), < tp < n, 
the imaginary part of Eq. can be rewritten as follows, 



cv 2 [2lnx + (2tp - 7r) cot tp] + 
+ [1 — vx{2 In >fCOS tp — (2tp — ir) sin tp)] 2 - 
+ [vx(2 In k sin tp + (2tp — ir) cos tp)] 2 = 0. 



(20) 



Starting from some threshold magnitude of x, there are 
two solutions of Eq. (|20|l for the phase tp at the given 
concentration of impurities, which correspond to the two 
existing bands. Respective values of e are then provided 
by the real part of Eq. lfT8)l. which closes up the paramet- 
ric solution of the problem. Correspondingly, the validity 
criterion for the CPA assumes the form, 
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Lax+l + i(p-f) 
lnx+ (tp — cot tp 



(21) 



As usual, for the renormalized wave vector in both 
bands one has ka = 2y/TTx\costp\. The spatial behav- 
ior of the host GF on one of sublattices at large intercell 
distances is given by 
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where Jo and Yq are the Bessel functions of the 1st and 
2nd kind, respectively. It follows from Eq. IL'l't that the 
mean free path should be written as I = a / (Ay/ir >c sin tp) . 
Thus, the localization parameter from the Ioffe-Regel 
criterionil takes the simple form, kl = | cot tp\/2. 

An overview of the SR scenario can be provided based 
on simple estimations. It may seem that e = a = cv 
is an appropriate solution of Eq. JTSJ. However, this is 
not the case. Formally, this equation is satisfied, but an 
analytical solution for the GF that is passing through 
this point can not be constructed. On the other hand, 
there should be always an energy, at which 5R(e — a) = 0. 
As follows from Eq. (|20|l . a certain amount of damping 
(5(7 7^ 0) is always present at this energy, 



2cv 2 In xq = 1 + (2v>cq In k^Y 



(23) 



When the impurity concentration is sufficiently small, 
kq w exp(— l/(2ct> 2 )). It is not difficult to see from the 
real part of Eq. (|18fl that 5Rcr w cv in this case. In other 
words, the energy, at which two bands coincide, is shifted 
approximately by cv from zero towards the impurity lo- 
cal level. The width of the concentration smearing area 
around e ~ cv, where states are highly localized accord- 
ing to the Ioffe-Regel criterion, should be proportional 
to kq, while the guess value for the mean free path in- 
side this area remains exponentially large. For small >r, 
Eq. J5D|) is reduced to 



ci> 2 (ln >c + (tt — 2tp) cot tp) — 1. 



(24) 



In the same approximation, it follows from the real part 
of Eq. (TH} that 



cv ~ cv 2 (tt — 2tp)>c/ sin tp. 



(25) 



Although the threshold magnitude of the localization pa- 
rameter, which separates states that can be described by 
the wave vector, can be argued to a known extent, it 
seems reasonable to choose it from the thoroughly tested 
method of potential-well analog y 12 ' 13 i 14 , |coty>| > v3. 
Then, the width of the concentration smearing area is 



*IR 



(8tt/3) exp(7r/v / 3)cz; 2 exp(-l/(2ci; 2 )). (26) 



In the narrow vicinity of e ~ cv, contributions from 
scatterings on impurity clusters are becoming signifi- 
cant. According to the applicability criterion Eq. I|21|l . 
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the electron spectrum obtained by the CPA can not 
be justified inside the area with the width of w 
exp(— l/(4ci; 2 ))/e, which is wider than Am that follows 
from the Ioffe-Regel criterion. It had been shown that in 
the 3D systems the small parameter of the series expan- 
sion Eq. I(17|) and the localization parameter kt can be 
expressed through each other 15 and depend on the phase 
(p only. However, in the system under consideration the 
cut-off phase for the CPA applicability criterion depends 
on the disorder parameter cv at the small impurity con- 
centration. The reason of this discrepancy is the subject 
of the more detailed study. 

With an increase in the impurity concentration, the 
absolute value of the shift \cv\ and the width of the con- 
centration smearing area are also gradually increasing 
in magnitude. It is obvious from the expression for xq, 
which can be rewritten as xq ~ — c|v|(2|i;|>fio lnxo) that 
parametrically \cv\ and xq simultaneously become of the 
order of the resonance energy |e r |. The second area of 
concentration smearing opens in the vicinity of the res- 
onance energy and, finally, both areas of concentration 
smearing are merged together. This is indicative of the 
spectrum rearrangement. Both criteria are coinciding in 
this regime (A# ~ A/a). An expression for the critical 
concentration of the SR can be obtained by comparing 
two main parameters of the problem by their magnitude, 

«; 2 exp[-l/(2CT 2 )] = (c\v\, (27) 

where £ is a certain constant to be determined. This 
immediately yields 

c r = -l/(2v 2 ln((/\v\)). (28) 

This expression fits well calculated critical concentrations 
of the SR with C~ KT 5 . 

At the impurity concentration that is far exceeding this 
critical value (i.e. c ^> c r ), it follows from Eq. 1)23(1 that 
in the first approximation xq does not depend on the 
magnitude of the impurity perturbation v . Both criteria 
give similar results for the width of the broad concentra- 
tion smearing area, A^ ~ Ajr ~ y/—(c)/ ln(i/c), which 
is nearly symmetric about the Dirac point of the host 



system. As was mentioned above, analogous approach to 
the description of the impurity effects in graphene have 
been undertaken in Ref. 0, in which some miscalcula- 
tions were committed in course of the theoretical analysis 
of the problem. Nevertheless, the width of the concen- 
tration smearing area for c 3> c r have been estimated 
correctly. 

While the passage of the SR in the system with linear 
dispersion deserves more close attention, some conclu- 
sions can be made at this stage. When the change in the 
on-site potential caused by the impurity atom is notice- 
ably larger than the bandwidth, the well-defined reso- 
nance state can do appear in the system with the Dirac 
spectrum. However, this resonance is not very sharp for 
the reasonable amount of the impurity perturbation. As 
a rule, the presence of the well-defined resonance state 
leads to the SR of the cross type with an increase in the 
impurity concentration. Yet, there are some exceptions 
from this rule, and the system under consideration be- 
longs to themi&. Despite the resonance, the SR is of the 
anomalous type that is common in low-dimensional sys- 
tems. This anomalous SR is characterized by the opening 
of the quasi-gap, in which any adequate cluster expan- 
sion can not be constructed and states are highly local- 
ized. The electron spectrum is not much distorted out- 
side of the concentration smearing area, and there are no 
prominent features in it close to the resonance energy. 

When the change in the on-site potential on the impu- 
rity site is not extremely large, the SR does not occur at 
all, and the width of the quasi-gap remains exponentially 
small. Indeed, from the practical point of view, such 
exponentially small quasi-gap will remain unnoticed in 
most situations, and virtually does not affect the carrier 
mobility. In the case of the large change in the on-cite po- 
tential disorder effects are not significant until the critical 
concentration of the impurities is reached. The obtained 
results also apply to the systems with the gap in the host 
quasi-particle spectrum when this gap is less than the 
width of the concentration smearing area. 

We are deeply grateful to Prof. V. P. Gusynin, who 
turned our attention to the problem of impurity states in 
systems with Dirac spectrum. 
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